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Abstract 

The scaling limit of the spin cluster boundaries of the Ising model with 
domain wall boundary conditions is SLE with k — 3. We hypothesise that 
the three-state Potts model with appropriate boundary conditions has 
spin cluster boundaries which are also SLE in the scaling limit, but with 
K = 10/3. To test this, we generate samples using the Wolff algorithm 
and test them against predictions of SLE: we examine the statistics of the 
Loewner driving function, estimate the fractal dimension and test against 
Schramm's formula. The results are in support of our hypothesis. 



1 Introduction 

Much of the recent interest in Schramm Loewner evolution (SLE) stems from 
the description it offers of cluster boundaries in critical statistical mechanics 
models. These models are conjecturally described by conformal field theories, 
which determines expectation values of products of local operators. However, 
non-local objects such as cluster boundaries do not have a natural expression 
in the field theory language. These curves are of great interest as, for example, 
they describe percolation cluster boundaries, level lines of height models and 
spin cluster boundaries. Recently, the scaling limit of the spin cluster boundary 
in the Ising model with appropriate boundary conditions has been proven to be 
SLE with K = 3 0. 

The Ising model corresponds to the g-state Potts model with q = 2. We may 
ask, therefore, whether spin boundaries in the critical Potts model with other 
values of q have SLE as their scaling limit. The Fortuin-Kastelyn (FK) cluster 
representation [?] of the Potts model contains closed loops which may be taken 
to be the contours of a suitably defined height model. The scaling limit of this 
model is supposed to be a field theory (CFT) called the Coulomb gas which 
possesses a single parameter g, where ^/q = — 2cos(7r(7). The connection with 
SLE is that the winding angle of these loops in the Coulomb gas are equivalent 
to that seen in radial SLE for k = 4/g. Therefore, we have reason to believe 
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that the FK cluster boundaries should correspond to SLE m the scalmg limit, 
with 

4 

K = -, r~n~ ■ (1) 

1 — arccos(y/^/2)/7r 

The spin cluster boundaries are related to these FK cluster boundaries in a non- 
trivial way, however, and it is therefore interesting to investigate the properties 
of these spin boundaries in the scaling limit. For the Ising model, these spin 
cluster boundaries are also described by SLE at the dual value n' = 16/k [S]. 

Besides — 1 (corresponding to percolation) , there are two other integer values 
of q for which the Potts model displays a continuous phase transition: q = i and 
(7 = 4 (for g > 4 the transition is first order). In this paper, we investigate the 
case g = 3 and test for agreement of numerically simulated lattice spin cluster 
boundaries with known predictions of SLE. The dual value of k for g = 3 is, 
from equation [1] above, k' = 16/k = 10/3. We may hypothesise, therefore, that 
the spin cluster boundaries, properly defined, of the q — i Potts model are SLE 
with K ~ 10/3. This paper is a numerical investigation of that prediction. 

The set-up of chordal SLE requires that the curves connect fixed points on the 
boundary of a simply connected domain. This is ensured by careful choice of the 
boundary conditions. For example, consider the following boundary conditions 
for the Ising model: divide the boundary into two connected subsets ('left' and 
'right'). Let the left boundary spins be spin up and the right boundary spins 
be spin down. This ensures that a spin cluster boundary runs from one fixed 
point on the boundary to another, for every configuration of the bulk spins. 
The Gibbs distribution induces a measure on these curves. The case of the 
g = 3 Potts model is similar, although there are now two types of interesting 
boundary conditions. If the boundary conditions are chosen as for the Ising 
model, with those on the left part being of spin type 1 (say) and those on the 
right of spin type 2, the right boundary curve of the cluster connected to the 
boundary spins of type 1 does not everywhere coincide with the right boundary 
curve of those of type 2. In other words, the spin cluster boundary may split 
around regions of spin type 3, which themselves may contain islands of the 
other spin types, ad infinitum. One possibility for the scaling limit is that the 
two cluster boundaries are distinct, with probability one. This might then be 
expected to correspond to two-curve SLE. We therefore simulate numerically a 
set of samples with these boundary conditions and compare to predictions from 
two-curve SLE: the fractal dimensions and the formulae describing the expected 
spatial distribution of the curves in the bulk [5] . The other choice for boundary 
conditions is to fix the left set of the boundary spins to be of type 1 and to allow 
those on the right to fluctuate between the other two spin types. In this case, 
we may consider the right boundary of the cluster connected to the boundary 
spins of type 1, which forms a curve connecting the two fixed points on the 
boundary. The measure on this set of curves is, as before, induced by the Gibbs 
distribution. Wc also simulate Potts models with these boundary conditions 
to test against two predictions of single-curve SLE: the fractal dimension and 
Schramm's formula. For the single curve case, we also examine the statistics 
of the driving function by undoing the sequence of conformal transformations 
which uniquely determine the form of the curve through the Loewner equation. 
This driving function should be Brownian motion with diffusivity k if the curves 
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correspond to SLE. 

The layout of the paper is as follows: In section [21 we summarise the definition 
and important features of the Potts model, explaining in detail the boundary 
conditions required to generate cluster boundaries. In section [31 we recall the 
definition of SLE and describe the predictions of SLE which we will test our 
samples against. In section[3]we describe the details of the numerical simulations 
we performed, explaining how the fractal dimension and other data are extracted 
from the generated samples. In section[51 we present the results of our numerical 
simulations and conclude in sectionlHl The appendix contains details of the Wolff 
algorithm used for all simulations. 

2 The Potts model 

The g-statc Potts model is a lattice spin model. The spins may take values 
from the set {1,2..., q} and have nearest neighbour interactions, such that the 
partition function takes the form 

Z = TTe'^ ^-y^-M-^'f-'K (2) 

J is the reduced coupling constant and 6ij is the Kronecker delta. The case 
5 = 2 is equivalent to the Ising model, up to a constant. For q > 2, therefore, the 
Potts model describes a generalisation of the Ising model to systems with more 
than two spin types. For a given g < 4 and choice of lattice, there exists a value 
of the coupling constant, Jc, at which the model has a continuous transition; the 
two point correlation function of the spins decays as a power law and the scaling 
limit of the model is supposed to be described by a conformal field theory. For 
later reference, these values for the two lattice types considered in this work 
are [U 

• Square lattice e'^" = y/q + 1. 

• Triangular lattice e"'^ = 2cos(| cos~^(y/g/2)). 

Let us consider spin cluster boundaries in the Ising model arising from a specific 
choice of boundary conditions. Namely, let those spins on a connected piece of 
the boundary be of spin type 1 and elsewhere on the boundary, the spins be of 
type 2. This is demonstrated in figure [1] where spin type 1 is represented by the 
colour red and type 2 spins are coloured green. Then, for every configuration 
of the spins in the bulk, there exists a unique curve starting and ending at two 
fixed points on the boundary. Immediately to one side of the curve, all the spins 
are green. To the other side, the spins are red. The measure on these curves is 
derived from the Gibbs distribution. 

Let us consider how to define the boundary conditions for the q ~ 3 Potts model 
in such a way as to guarantee the existence of cluster boundary curves propa- 
gating through the sample. The natural extension of the boundary conditions 
for the Ising model is to fix a connected set of the boundary spins to be of (say) 
spin type 1 (red or dark grey on black and white printouts), with the remaining 
spins of type 3 (blue or black) , as in figure [21 The particular choice of spins 
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Figure 1: The Ising model with boundary conditions as in the text. The spins 
are represented by hexagons, so that there are six nearest neighbours for each 
spin in the bulk. Half of the boundary spins are of type 1 (coloured red, or dark 
grey on black and white printouts) and the others of type 2 (green or light grey) . 
There exists a boundary of the cluster connected to the bottom edge (the thick 
black line) which coincides with the boundary of the cluster connected to the 
top edge. This curve connects the same two points on the boundary for every 
configuration of the bulk spins. 



is, of course, arbitrary due to the S3 symmetry of the model. We shall refer to 
these boundary conditions as 'fixed' boundary conditions. However, the situa- 
tion is different from the Ising model, since in this case the boundary curve of 
the cluster containing the red boundary spins docs not always coincide with the 
boundary of the cluster containing the blue boundary spins. In the figure, the 
dotted line is the boundary of the cluster containing the blue boundary spins, 
the dashed line is the boundary of the cluster containing the red boundary spins 
and where the two curves are both present along an edge, a thick black line has 
been drawn. This black line is reminiscent of the single cluster boundary in the 
Ising model, with spins immediately to one side of it red and those on the other 
side being blue, but the dotted and dashed lines are of a new type. The dotted 
line has blue spins to one side and either green (light grey on black and white 
printouts) or red to the other. The dashed line has red spins to one side and 
either green or blue to the other side. In this paper, we shall refer to the thick 
black lines as the 'composite' cluster boundaries and the dashed and dotted 
lines as the 'split' cluster boundaries. Note that the 'composite' cluster bound- 
ary splits into two 'split' cluster boundaries whenever it encounters a vertex 
with all three spin types, and that a pair of 'split' cluster boundaries recombine 
only at such three spin vertices. 

There exists another choice of boundary conditions for the q = 3 Potts model 
which we will refer to as 'fluctuating' boundary conditions. A connected set of 
the boundary spins are set to be of spin type 1 (red) . The remaining boundary 
spins are allowed to fluctuate with the spins in the bulk, but are not permitted 
to be of type 1. They are therefore of type 2 (green) or 3 (blue), see figure |31 
The cluster boundary of the red spin cluster connected to the top edge now 
coincides with the boundary of a cluster comprised of green and blue spins 
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Figure 2: The g = 3 Potts model 
with 'fixed' boundary conditions. 
The clusters connected to the bot- 
tom and top edges are no longer ev- 
erywhere overlapping. The dashed 
and dotted lines are 'split' cluster 
boundaries. The thick black line is 
the 'composite' cluster boundary. 



connected to the bottom edge, which is defined as follows: replace all green and 
blue spins with a new spin type x and consider the boundary of the cluster of x 
spins connected to the bottom edge. That these two cluster boundaries always 
coincide may be readily seen by analogy with configurations of the Ising model 
with spin types red and x. The reader will note the lack of reflection symmetry 
in the boundary conditions. This may be expected to lead to non- reflection 
symmetric spin cluster boundary curves. However, we shall present evidence 
that these curves arc actually reflection symmetric in the scaling limit. 



Figure 3: The q = i Potts model 
with 'fluctuating' boundary condi- 
tions. There exists a single uniquely 
defined boundary curve between the 
red (dark grey) cluster and the green 
(hght grey) and blue (black) cluster, 
as defined in the text. 



3 Schramm Loewner evolution 

We expect the cluster boundary curves of critical statistical models to be de- 
scribed by measures which are invariant under conformal transformations, due 
to the conformal invariance of the action associated to these models. This argu- 
ment can be made more precise; the changes in boundary conditions described 
in the previous section have been shown [21 [6] to be produced by the insertion 
of boundary 0i^2 operators in the appropriate conformal field theory. Expecta- 
tion values of observables therefore correspond to correlation functions with a 
pair of these boundary condition changing operators inserted. It may be shown, 
for a specific choice of observables, that the partial differential equations for 
these correlation functions arising from the null state associated with the 01^2 
operators are the same as those coming from the mathematical theory of ran- 
dom curves with conformally invariant measure, known as Schramm Loewner 
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evolution (SLE). 

Let us summarise the main features of the theoretical physicist's version OE] of 
SLE. A curve in the upper half plane from the origin to infinity is considered as 
being grown dynamically as a function of time t, which therefore parameterizes 
distance along the curve. Let Kt be the set consisting of the curve as grown up 
to time t (as well as all points enclosed by the curve and between the curve and 
the real axis) so that the complement of this set in the upper half plane is simply 
connected. Let gt{z) be the (unique) conformal mapping of this complement 
to the whole upper half plane, normalised such that gt{z) = z + 0{l/z) as 
z — > cxD. The coefficient of the 0(1/2) term is increasing with t, so 'time' can be 
reparametrised such that this coefficient is exactly 2t. Then, t is known as the 
capacity. The image of the growing tip of the curve under gt is a point at on the 
real axis. Loewner showed that the time-evolution of gt satisfies a differential 
equation, which is named after him: 

dgtjz) ^ 2 

dt gt{z) - at ' 

Any suitably continuous function at generates a curve. Schramm[7] showed 
that for this process to generate a conformally invariant measure on curves, the 
driving function must be a* = y/nBt + oq, with Bt being a standard Brownian 
motion. 

There remains a single parameter in the theory, k, which is the diffusivity of the 
Brownian motion. A small value of k results in a measure which favours curves 
with few turns, whereas a large value of k, leads to a measure favouring curves 
which wind more frequently, touching themselves and the real axis. The value 
of K is conjecturally related to the central charge of the conformal field theory 
it describes via the formula: 

_ (3^-8)(6-^) 

2k ■ ^' 

From this formula, it is clear that the central charge associated to an SLE with 
diffusivity k is equal to that of the dual SLE with k' = 16/k. 

We will now describe some properties of SLE curves. 



3.1 The fractal dimension 

If N{r) is the number of discs of radius r required to cover an object, this 
number scales as N{r) ~ r~'^f for small r, where the exponent d/ is known as 
the fractal dimension. SLE curves are fractal objects, with 

d/ = l + -^, for«;<8. (5) 
8 

The total length of the curve, S, measured in units of the lattice spacing, a, 
obeys a scaling relation as a function of system size, L, as L/a ^ 00 

S-af^y. (6) 
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3.2 Schramm's formula 



Schramm's formula relates to the following question: given an SLE connecting 
two boundary points, what is the probability that the curve passes to one side 
of a point, S,{z), in the interior? If the domain is mapped to the upper half 
plane such that the SLE connects to c», this probability may be expressed 
as the solution to a second order partial differential equation. Let us consider 
the probability that the curve passes to the left of a point, Picft- By this we 
mean that there is a continuous path from the point to the positive real axis 
which does not intersect the SLE curve. The solution is found by applying 
appropriate boundary conditions to the solution to the differential equation, 
which is a function only of t = Re(^)/Im(^) due to the conformal invariance of 
the measure: 



fle.w4-^^|^t2F,(i,l;|-i^). (7) 

The two domains considered in this work are related to the upper half plane by 
simple transformations, hence this result may be re-expressed easily in terms of 
the conventional coordinates in these domains. 

The generalisation of this result to the two-SLE case, was the subject of a 
recent paper [S]. Given a pair of SLEs both connecting to oo in the upper 
half plane, there are three possibilities: an interior point, S,{z), may be to the 
left of both curves, between them or to the right hand side of both. There is 
a probability for the curve to belong to each class, which is a function only of 
t = Re{^)/ Im{^), just as for the single curve case. The probability of each class 
obeys a third order linear partial differential equation, with solutions which 
arc integrals of hypcrgcomctric functions. The particular solution of interest is 
found by imposing a natural set of boundary conditions, namely the asymptotic 
behaviour as the point approaches the positive and negative real axis. The 
solutions are 



where 



and 



Pleltit) 



r(^)r(!) 



4.1. 



22-8/«7rr(i^ _ 1) 



2r(i+-4)r(A) 



r(i+10rFI+I) 



S{t')dt' , 



t2Fl{l 



i 3 
K ' 2 



4.3. 



(8) 



^'middlc(^) — 1 — 



r(^)r(!) 



22-s/K^Y{f - 1) 



Sit')dt' 



S{t')dt' 



(9) 



The probability that the curves are both to the right of a point may be deduced 
from the other two cases. Throughout this paper, we shall refer to these equa- 
tions as the two-curve formulae. As for the single curve formula, the results 
may be transferred to the domains of interest via a conformal transformation. 
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3.3 4-leg operators 



Consider the g = 3 Potts model with 'fixed' boundary conditions. There are a 
pair of cluster boundary curves in this case, which initially overlap, taking the 
form of a 'composite' cluster boundary. This may be expected to split around a 
cluster of spins which is not connected to the boundary. We may ask about the 
relevance, in a renormalization group sense, of the operator responsible for the 
recombination of these 'split' curves. This operator is familiar from the height 
representation of the 0{n) model, where it is known as the 4-lcg operator. The 
scaling dimension of this 4-leg operator is 

X4 = 2g- ^^-^ (10) 

where g = 1 — arccos (y/g/2)/7r. For q = 4, which we do not consider in this 
paper, the 4-leg operator is therefore marginal with X4 = 2, since 5 = 1 in 
this case. For q = 3, the scaling dimension is X4 = 99/60, which means it is 
irrelevant. As the system size increases, the 'split' boundary curves are therefore 
expected to collide less frequently. We may expect the scaling limit in this case 
to correspond to a pair of non-intersecting SLE curves. Note that the joining 
of two 'split' curves and a 'composite' curve at a point does not correspond to 
the insertion of a 3-leg operator, since the two types of curves are inequivalent. 
Indeed, we shall see that the 'composite' curve has a different fractal dimension. 



4 Numerical simulations 

We use the algorithm described in appendix |B] to generate a representative set 
of spin configurations at the critical point of the 5-state Potts model for both 
q = 2 and g = 3 in two domains. The first is the square lattice on the rectangle 
with aspect ratio 3:1, the second is the triangular lattice on the disc. 

As explained in section[21 if the boundary conditions for q ~ 2 (the Ising model) 
are chosen appropriately, a continuous boundary curve separating the two spin 
clusters runs between two points on the boundary. This ensemble of curves has 
recently been proven to converge to SLE3 [8]. It is therefore a useful check to 
show that the set of boundary curves obtained from our algorithm for q = 2 
satisfy the expected properties of SLE, as outlined in section [31 

The main aim of this paper is to gather numerical support for the spin bound- 
aries of the 3-state Potts model with appropriate boundary conditions becoming 
SLE in the scaling limit. There are two types of boundary conditions: 'fiuctuat- 
ing' boundary conditions (see section [5]) lead to a single curve running between 
two boundary points. 'Fixed' boundary conditions lead to a pair of 'split' curves 
running between two fixed points on the boundary, which may overlap to form 
a 'composite' curve. In section 1X51 we presented an argument for only the 'split' 
curves being present in the scaling limit. We test the validity of this argument 
by comparing the fractal dimensions of the two types of curves. 

The geometries chosen for the simulations are the square lattice on a narrow 
rectangular sample, and the triangular lattice on a disc. The first is chosen so 
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that Schramm's formula (see section [3. 2p may be compared with that on an in- 
finite strip. The second geometry is chosen for the ease of obtaining the driving 
function of single curves (see section 14. 4p . The two are complimentary in the 
sense that the measure on curves in the two domains should be related by the 
conformal transformation between them. That critical statistical models are 
conformally invariant is an accepted conjecture, but this does provide another 
useful check of the robustness of the algorithm, as well as a further demonstra- 
tion that properties of the scaling limit of the curves are consistent with those 
of SLE. 

The remainder of this section is dedicated to a technical discussion of the meth- 
ods of obtaining data from the samples. 

4.1 Identifying the cluster boundaries 

Given a configuration of spins, the first step to identifying the boundary curves 
is to identify the extent of the cluster containing a given boundary spin. This 
is done using the following algorithm: 

1. Start at the boundary spin of interest, marking it as belonging to the 
cluster 

2. Mark all neighbours if they are of the same spin type 

3. Repeat step 2 for all neighbours of marked spins 

4. Identify all connected sets of unmarked spins which are surrounded by a 
closed path of marked spins. Mark all such sets of spins. 

This process is demonstrated in figure U] for an Ising sample. 

It is sufficient to identify a single cluster boundary for the Ising model and q = 3 
model with 'fluctuating' boundary conditions, since the cluster boundaries of the 
top and bottom clusters overlap everywhere in these cases (see section [2|). For 
the case g = 3 with 'fluctuating' boundary conditions, as in figureOfor example, 
the first step to identifying the mixed spin cluster boundary is to replace all blue 
spins throughout the sample with green spins. Then, the algorithm is the same 
as that for the Ising model. 

For the <; = 3 model with 'fixed' boundary conditions, however, this process 
must be run twice to identify both the cluster connected to the bottom edge 
of the sample and that connected to the top edge, see figure \E\ The result 
is figure mf), which shows all three types of clusters present. Note that the 
example figures are for the triangular lattice case. The same algorithms to 
determine the cluster boundaries are employed for samples generated on the 
square lattice. 

4.2 Estimating the fractal dimension 

For the Ising model and the q ^ 3 Potts model with 'fluctuating' boundary 
conditions, there is a single curve running between two flxcd points on the 
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Figure 4: (a) shows a spin configuration in the Ising model. To obtain (b), 
recolour a red boundary spin at the top of the sample purple. Recolour all its 
red neighbours purple and repeat for all their neighbours et cetera. To obtain 

(c) , recolour a green spin on the bottom edge cyan. Recolour all its non-purple 
neighbours cyan and repeat for all their neighbours et cetera. Finally, to obtain 

(d) , recolour all remaining green and red spins purple. 

boundary. First, identify the cluster containing one set of boundary spins, 
as outlined in section 14.11 This results in a configuration of the type shown 
in figure UJd). Then, the length of the cluster boundary curve is deduced as 
follows. For each site on the sample, compare the colour at that site to each of its 
nearest neighbours. If the nearest neighbour is of the same colour (and therefore 
connected to the same cluster) there is not a cluster boundary separating the 
two spins. Otherwise, if it is of the other colour, add one half to the total length 
of the cluster boundary curve (since each link will be counted twice by this 
algorithm) . 

The situation for g = 3 with 'fixed' boundary conditions involves the fractal 
dimension of two types of cluster boundaries: the 'composite' cluster boundaries 
and the 'split' cluster boimdaries, see section[51 Identify the clusters containing 
the two types of boundary spin and also the clusters which do not contain a 
boundary spin. This leads to a cluster configuration similar to that shown in 
figure[51Jf). The total length of 'composite' boundaries is found by summing over 
all sites, adding one half to a running total wherever a pink site is adjacent to a 
grey site. The total length of 'split' cluster boundaries is found by identifying all 
sites where a grey site is adjacent to a yellow site or a pink site is adjacent to a 
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Figure 5: Figure (a) shows a g = 3 Potts model spin configuration with 'fixed' 
boundary conditions. To obtain (b), recolour a red boundary spin at the top 
of the sample purple. Recolour all its red neighbours purple and their red 
neighbours purple et cetera. To obtain (c), recolour a blue spin on the bottom 
edge cyan and all its non-purple neighbours also. Repeat for all their non- 
purple neighbours et cetera. Recolour purple all the remaining red, blue and 
green spins. To identify the cluster containing the blue boundary spins, start 
again from figure (a). Fig (d) shows the result of starting with a blue spin on 
the bottom edge and recolouring all its blue neighbours purple, et cetera. Fig 
(e) is obtained after recolouring cyan all non-purple neighbours starting from a 
red spin at the top edge. Lastly, recolour all remaining red, blue and green spins 
purple. The two clusters connected to the boundary are coloured differently in 
figures (c) and (e). But some spins are cyan in both (c) and (e). These are 
yellow in figure (f). The top cluster is pink and the bottom cluster is grey. 

yellow site, adding one half to a different running total for each such occurence. 

In this way, we may deduce the expected length of curves as a function of system 
size. Comparison to equation [6] allows an estimation of the fractal dimension of 
the curves. 

4.3 Comparison to Schramm's formula and its generalisa- 
tion to two curves 

Firstly, for the Ising model and q = 3 Potts model with 'fluctuating' bound- 
ary conditions, the starting point is the cluster configuration picture, as in 
figure m^d). Recall that the boundary conditions are chosen such that the curve 
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starts and ends at opposite points on the boundary. There is a path which is 
equidistant from these two points, which is shown in figure [S] for the rectangle. 
For each site along this path, calculate the percentage of the samples for which 
the site is coloured cyan. This yields the fraction of samples for which this point 
is connected to the bottom cluster. This may be compared to the predictions 
of Schramm's formula. 

For the q = 3 Potts model with 'fixed' boundary conditions, the procedure is 
to deduce the cluster picture, as in figure [SUf). Two separate quantities may be 
measured, the percentage of the samples for which each point along the path is 
connected to the bottom cluster (coloured grey) and the percentage which are 
in the central, yellow cluster. In the former case, the point is to one side of both 
curves (the left hand side say), whereas in the latter example, the points are 
between the two curves. These results are to be compared to the generalisation 
of Schramm's formula to two curves, the two-curve formulae (see section 132]) . 



4.4 The driving function 

The third test is to analyse the driving function of the curves for the Ising model 
and q = 3 Potts model samples with 'fiuctuating' boundary conditions. The spin 
boundary curves have a natural discretisation provided by the lattice itself, so a 
curve may be approximated by the set of ordered lattice points which it passes 
through. It is easiest to consider the set of conformal mappings of the curve 
back to the boundary in the upper half plane, so the first step is to transform 
domains using a conformal transformation. For the disc, for example, this is 
the following Mobius transformation: 

For the rectangle, the transformation to the upper half plane is not so conve- 
nient; the derivative of the conformal map at the corners is large. The effect of 
this cannot be simply and systematically removed. 

The curve grown up to a point 7 may be mapped back to the real axis by a 
series of vertical slit maps for each successive discretizcd point. Xii along the 
curve 

g,,t(z) = Re(xO + Vii^ - Re(xO)2 -f Im(x.)2) . (12) 

This maps the point Xi to Re(xi) on the real axis, which is therefore the value 
of the driving function, at, at the time when this point is the tip of the growing 
curve. The change in t associated with this piece of curve is found from the 
expansion of the map as ^ 00 : 

Im(x)V2 
gt{z ^ 00) = z + — ' + ... 

z 

The change in t is therefore St ~ Im(xi)^/4. For SLE, the driving function 
should be a multiple of Brownian motion. In particular, the expectation value 
of at should vanish and the variance should be k. On a plot of frequency of step 
lengths (normalised by ^/St) the distribution should be Gaussian with variance 
Kt. Our lattice curves are only supposed to approximate SLE at scales much 
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larger than the lattice spacing. Therefore, effective steps equivalent to 10 lattice 
steps are considered. SB is then the change in the driving function as a result of 
an effective step, and 6t is the capacity change, totalled over the 10 component 
steps. We also plot (a^ ) as a function of i, with a view to equating the gradient 
to K. 



5 Results 

In this section, we present a comparison between the results from the numerical 
simulations and predictions from SLE. The algorithm used to generate samples 
is described in appendix [Bl 

5.1 Ising model on the rectangle 

This section describes the results from the Ising model on the rectangle on a 
square lattice of aspect ratio 1:3. Four system sizes are used: 40 x 120, 60 x 180, 
80 X 240 and 100 x 300, where the lengths are in multiples of the lattice spacing. 
These will be referred to by their smaller dimension. The boundary conditions 
arc as described in section [21 with one spin type fixed along a long edge and 
half each of the shorter sides, and the other spin type elsewhere. The geometry 
is therefore as shown in figure El A cluster boundary curve runs from the 
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Figure 6: The geometry used for the Ising model on the rectangle of aspect ratio 
1 : 3 with a square lattice. The boundary conditions are shown, but the bulk 
spins left grey. The path used to test against Schramm's formula is shown as 
the blue dashed line. 

centre of one short edge of the rectangle to the centre of the other short edge. 
50, 000 independent samples of each system size are used to estimate the fractal 
dimension of the curves, by determining the average length as a function of 
system size. This process is explained in section 14.21 The results are plotted 
in figure [71 The error bars are shown, but are only just discernible at this 
magnification. The linear fit is in good agreement with the data, with reduced 
equal to 0.7. It shows a fractal dimension of df = 1.372 ± 0.002. Using 
formula El this would correspond to SLE with k = 2.976 ± 0.02. The known 
value of K = 3 for the Ising model lies just outside these error bars. 

The fit to Schramm's formula is via the logarithmic map from the upper half 
plane to the infinite strip: 

, , In z , ^ 

w{z)^—. (13) 
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Fractal dimension of the spin oluster boundary curve tor ttie Ising model on the rectangle 



numehcal data ^ 
y.1.372'x+0.717 



1.6 1.65 



1.75 1.8 1.85 

ln(system size} 



1.9 1.95 



Figure 7: The fractal dimension of the Ising model on the rectangle. The system 
size is defined as the number of lattice spacings across the smaller side of the 
rectangle. The curve length is measured in multiples of lattice spacings. 



The assumption that the strip is effectively infinitely long is valid only if the 
boundary curve at the centre of the sample is independent of the boundary 
conditions at the ends of the rectangle. We tested this by generating system 
sizes of aspect ratios 1 : 5 and 1:6. The differences in the data are smaller than 
the error in the measurements. In terms of a coordinate, s, which describes 
the fractional distance across the strip, t is given by t = cot(7rs) in equation [71 
This follows from equation [T31 The results from the numerical data are the 
expectation value for the point at a fractional distance s across the strip being 
connected to the bottom spin cluster. Figure [8] shows the comparison to the best 
fit value of k = 3.021 ± 0.015. Once more, the known value lies just outside this 
range. We should not be surprised, therefore, it the predicted value of 10/3 for 
9 = 3 Potts model also lies just outside the range obtained from the numerics. 
The average is again over 50,000 independent samples. 



5.2 g = 3 on the rectangle 

The geometry is the same as for the Ising model above, with sample sizes 40 x 
120, 60x 180, 80x240 and 100x300. First, let us examine the case of 'fluctuating' 
boundary conditions, when there is a single curve connecting two boundary 
points at the centre of each short end of the rectangle. The boundary conditions 
are like those in figure El except that the boundary spins on the bottom half are 
permitted to be either blue or green and to fluctuate with the spins in the 
bulk. The fractal dimension comes from the average length of the curves as a 
function of system size, where 50, 000 independent samples have been generated 
for each sample size. The results are shown in figure [5] The linear fit shows 
the fractal dimension to be d/ ~ 1.399 ± 0.002, which would correspond to SLE 
with K ~ 3.192 ± 0.02. The fit to Schramm's formula, figure \10\ shows that 
on average, the curve is closer to the fluctuating spins boundary than would 
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Agreement with Sctiramms tormula for ttie Ising model on the rectangle 




0.2 0.4 0.6 0.8 1 

Distance of point across the sampie 

Figure 8: Fit to Schramm's formula on the strip for the Ising model, system size 
100 X 300. The smaller system sizes also show good agreement with Schramm's 
formula with k w 3.0. The agreement is to approximately 1%, as shown in the 
lower figure. 



be expected. Equivalently, a point of fractional distance < s < 1 across 
the sample from the fluctuating boundary towards the fixed boundary is more 
likely to be connected to the fixed spin cluster than would be expected from 
Schramm's formula. Perhaps this is not surprising; whereas the two types of 
boundary conditions in the Ising model are reflection symmetric, the boundary 
conditions for g = 3 with fluctuating boundary conditions are not. It is only 
in the scaling limit that we might expect to recover reflection symmetry. The 
best fit value is k = 3.235 ± 0.01. We shall see that this asymmetry is a result 
of non-zero average of the driving function, but that this effect is localised to 
the neighbourhood of the boundary spins, and therefore is unimportant in the 
scaling limit. 

An alternative set of boundary conditions are the 'fluctuating' boundary condi- 
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Fractal dimension for q=3 Potts modei witti 'fiuctuating' boundary conditions on the rectangie 
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Figure 9: The fractal dimension of the q = 3 Potts model on the rectangle with 
fiuctuating boundary conditions. System size is the shorter side of the rectangle 
in units of lattice spacing. The curve length is in units of lattice spacings. 



tions described in section[21 These lead to the existence of two cluster boundary 
curves; the boundary of the cluster containing one set of boundary spins and 
the boundary of the cluster containing the other boundary spins. Wherever the 
two curves overlap, this is labelled a 'composite' part of the curve. Where they 
are distinct they arc labelled 'split' curves. In order to ascertain the fractal di- 
mensions of the two types of curves, the expected total length of the two types 
of curves as a function of system size is plotted, see figure [TT] 65, 000 indepen- 
dent samples are generated for each sample size. The fractal dimension of the 
'composite' boundary type is df = 1.022 ± 0.003, which would correspond to 
SLE with li = 0.18±0.02. The reduced of this linear fit is 9.7 however, which 
is indicative of non-linearity in the results. Closer investigation reveals convex- 
ity resulting from the fractal dimension diminishing as system size increases. 
The 'split' boundary type has a larger fractal dimension of df = 1.598 ± 0.008, 
which would result from SLE with k = 4.78 ± 0.06. This supports the hypoth- 
esis that the 'split' boundary types dominate in the scaling limit, resulting in a 
pair of curves running between the boundary points. The value of k is a little 
surprising; it is close to the dual value for k = 10/3, namely k' = 16/k = 4.8. 

Figure [T^l shows the agreement with the generalisation of Schramm's formula 
to two curves, again from 65, 000 independent samples. The path used is the 
same as for the single curve case, namely the dashed line in figure [6l As a 
function of the distance of a point upwards from the bottom boundary, we may 
determine the expected fraction of the samples for which both cluster boundary 
curves arc above the point (top figure) and the fraction for which the point 
is between the two curves (bottom figure). The fit to the two-curve formula 
is poor for the lattice sizes considered, as may have expected; for finite lattice 
size, there is a non-zero contribution from 'composite' cluster boundaries, which 
reduce the proportion of points which lie between the two curves and increase 
the proportion of points to one side of both curves as a result. 



16 



Schramm's formula for q=3 Potts model on file rectangle witfi 'fluctuating' boundary conditions 



numerical data ^ 
kappa=3.235 



Distance of point across tfie sample from fluctuating boundary 
Difference between numerical data and Schramm's formula 




Distance of point across the sample from fluctuating boundary 



Figure 10: Schramm's formula for q = 3 on the rectangle for system size 100 x 
300. The distance across the sample is from the fluctuating boundary toward 
the spin 1 boundary. 



5.3 Ising model on the disc 

The spins lie on the triangular lattice so that each spin has six nearest neighbours 
and the spin boundaries lie on the edges of the honeycomb lattice (the sides of 
the hexagons). Initially, four system sizes are considered, characterised by the 
maximum number of rows across the central vertical diameter: 101, 153, 201 
and 253. Spins are included if the centre of the hexagon is a distance away 
from the centre of less than or equal to the radius of this vertical column. The 
geometry is shown in figure [T3l with the boundary conditions used for the Ising 
model. The fractal dimension of the curve in this geometry is df = 1.38± 0.002, 
which would correspond to SLE with k = 3.08 ± 0.02, see figure [HI This result 
is obtained from 10, 000 independent samples for each system size. 

The fit to Schramm's formula is via the Mobius transformation from the upper 
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Fractal dimension of the q=3 Potts model 'composite' boundary on the rectangle 



numerical data ^■ 
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Fractal dimension of q=3 Potts model 'split' boundary on the rectangle 



numerical data ^■ 
1 .598'x+0.425 
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Figure 11: Fractal dimension of 'composite' and 'split' boundary types for q — 3 
Potts model on the rectangle with fixed boundary conditions. The system size 
is the length of the smaller side of the rectangle in units of lattice spacing. 



half plane to the disc 

= (14) 

The value of t in equation [7] is related to the fractional distance < s < 1 
upwards along the vertical diameter by the formula 

2s - 1 
2s(l — s) 

The best fit for the largest system size considered would correspond to SLE with 
K = 3.018±0.007, as shown in figure [T5l The reason that only the largest system 
size is shown is that the plot is quantitatively similar for all system sizes. Over 
the range of system sizes tested, the value of k, fluctuates around a constant 
value. This appears inconsistent with the hypothesis that the value should be 
converging to k = 3 as system size increases, but is likely to be a result of 
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Comparison to two-curve formula for q=3 Potts model on the rectangle 
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Figure 12: Comparison to the two-curve formula for the three-state Potts model 
on the rectangle. The agreement with the prediction oi k ~ 4.8 from the fractal 
dimension is poor. Two system sizes are shown, / = 40 and / = 100, where / is 
the smaller of the sides of the rectangle. 



the small range of system sizes considered. Larger systems are computationally 
demanding to simulate, but would offer further insight into this problem. 

The third test against SLE is to extract driving functions for each curve and to 
examine their statistics as a function of system size. This process is described in 
section An example, for diameter 253, is shown in figure [TOl The results are 
from 50, 000 independent samples. For each sample, only the part of the curve 
corresponding to capacity t < 1/3 was used, since the error introduced by the 
series of conformal maps may be expected to grow exponentially. This would 
result in inaccurate values for the driving function from later points on the curve. 
The form should be Gaussian, but the finite lattice is apparently responsible 
for a bi-modality for small steps SBt/VSi. A second test is to examine the 
expected squared value of the driving function as a function of capacity, t. Each 
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Figure 13: The geometry for the Ising model on the disc. The figure would 
be classified as size 13 as this is the number of spins across the largest vertical 
diameter. Boundary spins are fixed as shown and bulk spins left grey. The path 
used for the comparison to Schramm's formula is shown as the dashed cyan line. 



Fractal dimension of the ising modei on the disk witii a triangular lattice 



2.9 - 



2.8 I ' 1 1 ' ' ' 1 ' 1 

2 2.05 2.1 2.15 2.2 2.25 2 3 2 35 2.4 2.45 

ln(system size} 

Figure 14: The fractal dimension of the Ising model on the disc with a triangular 
lattice from 10, 000 samples of each system size. The size of the system is the 
number of spins across the central vertical diameter. 



system size has a linear fit with a gradient which would correspond to k. The 
results for diameter 253 is shown in figure 1171 The linear fit is excellent, with 
reduced = 0.01. The value of this best fit as a function of the inverse system 
size is plotted in figure 1181 The approach to the scaling limit appears to be 
complicated; it is not a simple power law at these lattice sizes. 

Smaller system sizes and larger numbers of independent samples are used for 
this test, as large numbers of samples are required to obtain results to a desirable 
accuracy. The table below shows the system sizes used, with the corresponding 
values of the inverse number of spins and number of samples generated. 
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Comparison to Schramm's formula for Ising model on the disk diameter 101 




-0.004 I 1 ^ ^ 1 1 

0.2 0.4 0.6 0.8 1 
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Figure 15: The best fit to Schramm's formula for the Ising model on the disc 
with 101 spins across the vertical diameter.. The agreement is to better than 
1%, see the lower figure. 



System size 


Inverse number of spins 


Number of samples used 


101 


1.47 X 10-4 


760,000 


113 


1.18 X IQ-^ 


920,000 


125 


9.58 X 10-5 


720,000 


137 


7.98 X lO-'^ 


1,040,000 


153 


6.37 X 10-5 


580,000 



5.4 q = 3 Potts model on the disc 

Firstly, we examine the case of 'fluctuating' boundary conditions, such that 
a single cluster boundary propagates through each sample. System sizes are 
again characterised by the number of spins along the central vertical diameter. 
System sizes considered are 101, 153, 201 and 253. The fractal dimension. 
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step distribution from the driving function for the Ising model on disk radius 253. 
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Figure 16: The frequency of step distribution for the Ising model on the disc, 
diameter 253. The plot takes data from 50,000 independent samples. The 
y-scale is an arbitrary measure of frequency. 



obtained as usual by determining the expected curve length as a function of 
system size, is plotted in figure fT^ The value for the fractal dimension obtained 
is df = 1.4216 ±0.0007, which would correspond to SLE with k = 3.373 ±0.006. 
This comes from considering 10, 000 independent samples of each system size. 
The fit to Schramm's formula, as for the Ising model, fluctuates with system 
size across the range of sample sizes considered. For the largest system size 
considered, the best fit is to the predictions from SLE with k = 3.275 ± 0.003. 
This is shown in figure !^ The driving function may be obtained by considering 
the series of conformal maps which undo the curve up to capacity, t, as for the 
Ising model in the previous subsection. A similar agreement with a Gaussian fit 
is again seen, with the same characteristic bi- modal behaviour at small SBt/\/Si. 
It is more instructive to examine (B^) as a function of t (the variance of the 
distribution) and to extract the gradient of the linear fit. Each system size shows 
linear behaviour to a high degree of accuracy. A typical plot is as in figure [2T| 
which has reduced = 0.01. The values of the gradient of the linear fit to 
(flj) against capacity, as a function of inverse system size, is shown in figure !^ 
As for the Ising model, smaller system sizes arc used for this test because large 
number of independent samples are required to keep the error manageable. As 
for the Ising model, the system does not show a power law approach to the 
scaling limit at these system sizes. 



System size 


Inverse number of spins 


Number of samples used 


101 


1.47 X 10-^ 


800,000 


113 


1.18 X 10"^ 


800,000 


125 


9.58 X 10-5 


800,000 


137 


7.98 X 10-5 


640,000 


153 


6.37 X 10-5 


720,000 



Recall that, by comparison to the fit to Schramm's formula, the spin boundary 
curve of the g = 3 Potts model with 'fluctuating' boundary conditions displays 
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Mean square of the driving function for the Ising model on disc diameter 253 



numerical data ^■ 
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Figure 17: The expected value of plotted against capacity, t, for the Ising 
model on the disc diameter 253. The fit is excellent, with reduced = 0.01. 
The data comes from 40,000 independent samples. 



asymmetry. The curve is more likely to pass close to the fluctuating boundary 
than is predicted by SLE. This asymmetry is not surprising; the boundary 
conditions are not reflection symmetric, as they are in the Ising model. However, 
this asymmetry must vanish in the scaling limit if the curve is to become SLE. 
Figure [23] shows the mean of the driving function for the curve on the disc of 
diameter 253 as a function of capacity, t. It is clear that the asymmetry is 
a result of the behaviour of the curve at small t. This is encouraging, since 
it suggests that this is a boundary effect, whose contribution would therefore 
become irrelevant in the scaling limit. 

The case of 'fixed' boundary conditions may be tested for fractal dimension of 
the two types of boundary, the 'composite' and 'split' boundaries. A plot of the 
fractal dimension of each of these is shown in figure [211 40, 000 independent 
samples of each system size were used. The fractal dimension of 'composite' 
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Figure 18: The mean gradient of the hnear fit of to t for the Ising model as 
a function of system size. The fit shown is hnear, which is clearly inadequate; 
the approach to the scaling limit is not a simple power law at these scales. 



Fractal dimension of q=3 Potts l\^odel on disk radius 253 
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Figure 19: The fractal dimension of the q = 3 Potts model on the disc. System 
size is the number of spins across the central vertical diameter. 



curves is not constant but diminishing as system size increases, just as was 
seen for the model on the rectangle. The fractal dimension of the 'split' curves 
df ^ 1.589 ± 0.005 would correspond to SLE with k = 4.71 ± 0.04. The fractal 
dimension of the the 'composite' curves is df = 1.03 ± 0.01 {k — 0.24 ± 0.08). 
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Comparison to Schramm's crossing formuia for q-3 Pctts modei on the disk. 
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Figure 20: Best fit to Schramm's formula for the g = 3 Potts model on the disc 
diameter 101. The same asymmetry is seen as for the model on the rectangle. 



6 Summary and Conclusions 

Recall our hypothesis, that the scaling limit of spin cluster boundaries in the 
three-state Potts model is SLE with k — 10/3. The numerical results from the 
simulations with 'fluctuating' boundary conditions, and hence with single curves 
propagating across the samples, shows encouraging support for this hypothesis. 
The data should be compared to that from the Ising model, which has been 
proven to correspond to SLE^. The results are summarised in the table below 
for the tests of the fractal dimension and best fit to Schramm's formula, both 
on the rectangle and the disc. 



Rectangle 


Ising 


g = 3 


Fractal dimension 
Schramm's formula 


2.976 ±0.02 
3.02 ±0.02 


3.192 ±0.02 
3.235 ±0.01 
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Mean square of the driving function for the q=3 Potts model on the disc 
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Figure 21: The mean square value of the driving function as a function of 
capacity, t. The system is the disc of diameter 201 and 40, 000 independent 
samples are used. Agreement with the linear fit is excellent. 



disc 


Ising 


q = 3 


Fractal dimension 
Schramm's formula 


3.08 ±0.02 
3.018 ±0.007 


3.373 ±0.006 
3.275 ±0.003 



The value for the fractal dimension is the more robust result, since it incorpo- 
rates data from the whole range of system sizes. Notice that the known value 
of K = 3 lies just outside the range obtained from the fractal dimension of the 
numerical simulations on the rectangle. The values obtained are a little too 
low for the Ising model, but also for the three-state Potts model, compared to 
our prediction oi k = 10/3. However, on the disc, the Ising model shows a 
fractal dimension which is larger than the known value. In this case, the fractal 
dimension of the Potts model is also larger than our prediction. Qualitative 
agreement with Schramm's formula for a single system size (100 x 300 on the 
rectangle and diameter 253 on the disc) is also encouraging. We have demon- 
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Figure 22: The gradient of the linear fit to the mean square value of the driving 
function against capacity, plotted against inverse system size. The approach to 
the scaling limit is not a simple power law. 



IVIean of ffie driving function against capacity for tfie q=3 Potts model on tfie disc 




0.04 0.06 
capacity, t 



Figure 23: Mean of the driving function as a function of capacity for the g = 3 
Potts model on the disc diameter 253. 



strated that the asymmetry of the curve (it is more likely to be close to the 
fluctuating boundary than the prediction from SLE) is a result of the behaviour 
of the curve at small capacity, as may be seen from a plot of the mean of the 
driving function against capacity. This asymmetry may be assumed, therefore, 
to be a boundary effect, which would be irrelevant in the scaling limit. This 
infers that reflection symmetry is recovered in the scaling limit, consistent with 
the curves being SLE. 

Since the results from the variance of the driving functions do not appear to 
approach the scaling limit as a simple power law, it is not clear how to extract 
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Fractal dimension of composite boundat7 cun/es for q=3 Potts modei on tfie Disk 
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Figure 24: The fractal dimensions of the 'composite' (top figure) and 'spUt' 
(bottom figure) boundaries for the q = 3 Potts model on the disc. 



a value for k from the data. If we simply average over the values obtained from 
all system sizes considered, we find k ~ 3.04 for the Ising model and k ~ 3.30 
for the g = 3 Potts model. This figure is, of course, an arbitrary one and is 
supposed only to give a guide to the order of magnitude of the variance in the 
scaling limit. The fact that this figure is close to the predicted value is, however, 
encouraging. 

We also considered the three-state Potts model with fixed boundary conditions, 
which leads to a pair of cluster boundaries, coinciding in some places. The 
results for k from the fractal dimension of these two curve types on the rectangle 
and the disc are summarised below. 



2 Curves 


Composite 


Spht 


Rectangle 
disc 


0.18 ±0.02 
0.24 ±0.08 


4.78 ±0.06 
4.71 ±0.04 
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The 'split' boundaries have a larger fractal dimension. This is consistent with 
the hypothesis that they dominate in the scaling limit, supported by the irrel- 
evance of the 4-leg operator. For the sample sizes considered, the agreement 
with the two-curve formula (the generalisation of Schramm's formula to two- 
curve SLE) is poor. This is unsurprising, as the presence of 'composite' cluster 
boundaries reduces the frequency of points being between the two curves and 
redistributes this to the probability that a point is to one side of 'both' curves. 
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A The Fortuin-Kastelyn cluster representation 



There exists a well known mapping of the partition function to a cluster model, 
which will be referred to in the appendix describing the Wolff algorithm. Up to 
a multiplicative factor, the partition function is 

Z = TtY[{^ ~ P + P6s(r),s(r')) ' (16) 
r,r' 

with e~^'^ = 1 — p. The product may be expanded into a sum of 2^ terms (N is 
the number of nearest neighbours) by choosing either (1 — p) or pSs(^r),s{r') from 
each bracket. If a line is drawn on the edge between each pair of spins for which 
the term p(5s(r),s(r') "^^s chosen, the partition function is seen to correspond to 
a sum over all graphs on the edges of the original lattice. The spins in each 
connected cluster are constrained to take the same value, so taking the trace 
leads to _ 

^ = ^p|G|(i_p)|G|g||G||_ (17) 
G 

The sum is over all graphs, \G\ is the number of edges in a given graph, \G\ is the 
number in the complement and \\G\\ is the number of connected components, 
known as Fortuin-Kastelyn (FK) clusters. 



B The Wolff algorithm 



The Monte Carlo algorithm used throughout to generate samples was the Wolff 
algorithm fTUl E]. It is a Monte Carlo algorithm with the following update 
procedure: 

1. Pick a lattice site at random. 

2. Visit all neighbouring sites, adding them to the cluster with a probability 
p if the spin is of the same type. 

3. Repeat step (2) for all the newly included sites, adding their neighbours 
with the same probability if they are also of the same type. 

4. Repeat (3) until no new sites may be added. 

5. Change all spins in the cluster to a randomly chosen new spin type. 



There is always a sequence of single spin flips from any configuration to any 
other. The probability, p, is related to the reduced coupling by (cf equation III 



e^-^ - 1 

It is the conditional probability that two neighbouring spins of the same type 
belong to the same FK cluster. The algorithm is therefore equivalent to iso- 
lating a single FK cluster and changing all included spins to a new type. The 
benefit of using this algorithm is that the phenomenon of critical slowing down 
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is eliminated. However, the average computational time required for an update 
increases as the critical temperature is approached due to increasing average 
cluster sizes. 

This procedure describes a single update. To generate a sample, we start with 
a random initial configurations of spins which are consistent with the boundary 
conditions. Then the update algorithm is run a sufficient number of times for the 
energy of the system to equilibriate. In this way, we arrive at the first sample. 
To calculate the number of updates required to produce a second, independent 
sample, we examine the autocorrelation function: 

u i 

for a number of samples such that the average over ti yields convergent results 
for large enough tj. Note that (si) is also a time average. This autocorrelation 
function decays exponentially as G{tj) ~ e~*^^'^. The correlation time, r, is a 
function of system size. Samples are deemed independent after 3 correlation 
times. 

In this paper, we consider two boundary types. The first is fixed to spin type 
1 (say). If a cluster generated by the Wolff algorithm connects with such a 
boundary, we do not flip the spins. That is to say we ignore step (5), thus 
keeping the boundary spins fixed to the only allowed type. The other choice of 
boundary condition considered is for the spins not to be allowed to be spin type 1 
(say) . If a cluster generated by the algorithm connects with a boundary of this 
type, then the allowed spin flips in step (5) of the algorithm are constrained to 
not be spin type 1. In this way, the boundary spins fluctuate with the sample, 
but never take the value 1. 
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